Discreteness and entropic fluctuations in GREM-like systems 
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Within generalized random energy models, we study the effects of energy discreteness and of 
entropy extensivity in the low temperature phase. At zero temperature, discreteness of the energy 
induces replica symmetry breaking, in contrast to the continuous case where the ground state is 
unique. However, when the ground state energy has an extensive entropy, the distribution of overlaps 
P(q) instead tends towards a single delta function in the large volume limit. Considering now the 
whole frozen phase, we find that P(q) varies continuously with temperature, and that state-to-state 
fluctuations of entropy wash out the differences between the discrete and continuous energy models. 



I. INTRODUCTION 

Strongly disordered systems such as spin glasses^ have 
many metastable states in their low temperature frozen 
phase, rendering them very sensitive to perturbations. 
An extreme case of sensitivity arises in systems having 
replica symmetry breaking) 2 ! (RSB): there the excess free- 
energy of the lowest excited states is O(l) and so P(q), 
the distribution of overlaps q, is non-trivial. Given such 
sensitivity, what properties are robust to changing the 
details of the microscopic Hamiltonian, and furthermore 
can the presence of RSB itself depend on microscopic 
details? Our purpose is to investigate this point in the 
context of tractable models of spin glasses which are of 
the mean field type. The main motivation for this is the 
question of RSB at T = in the ± J Ising spin glass. It 
has been argued by Krzakala and Martina that in phys- 
ical systems with highly degenerate ground states there 
should be no replica symmetry breaking at T = even if 
there is RSB at T > 0. Numerical investigations of this 
issue in the three-dimensional Ising spin glass have led to 
conflicting claims, either validatingErU this picture or on 
the contrary^ suggesting that there is RSB amongst the 
ground states of that model- To understand better this 
question, Drossel and MooreO have investigated overlaps 
in the ± J model on the Migdal-Kadanoff lattice at T = 
and T > 0, concluding that P{q — 0) behaves differently 
with lattice size in the two cases; however that model 
does not have RSB at any temperature. In this work 
we provide a study of this question in a solvable model 
having RSB in its low temperature phase. n 

This paper begins with the random energy modelEl of 
Derrida and we consider the effects of having discrete en- 
ergies. Because the ground state can be degenerate, P(q) 
at T = has a strictly positive weight at q — while in 
the continuous case P(q = 0) goes to zero linearly with 
temperature- Then we consider generalized random en- 
ergy modelsu and in particular a discrete version with an 
infinite number of layers that can be compared tp-the con- 
tinuous version studied by Derrida and Spohriid. In all 



these cases, the discreteness gives rise to RSB at T = 0. 
One of the unphysical features of such models is their 
zero entropy density at low enough temperature. Since 
we expect entropy fluctuations to be important for the 
question of replica symmetry breaking, we extend this 
last model so that states have random entropies. For 
the models considered with an infinite number of layers, 
we compute P(q), the disorder averaged probability dis- 
tribution of overlaps. We then see the effects of energy 
discreteness and of entropy fluctuations, in particular as 
T — > 0. Finally, we consider how replica symmetry is 
restored at T — in the thermodynamic limit through 
entropic fluctuations, even though the energies are dis- 
crete. 



II. RANDOM ENERGY MODEL 

The REm! is a simple model of spin glasses; its parti- 
tion function is 

^=i>p[-§] a) 

i=l 

where N is identified with the number of spins of the 
system and the units are chosen so that the Boltzmann 
constant is 1. The energies -E 1 , are independent iden- 
tically distributed random variables of law Pn{E). To 
make contact with physical systems, Pn(E) is set to co- 
incide with the distribution of energies of spin configu- 
rations in the d-dimensional Edwards-Andersonliil (EA) 
spin glass. Such a system with N spins has 2 N energy 
levels but these energies are strongly correlated; the REM 
is obtained if one neglects these correlations. If the cou- 
plings Jij of the EA model are Gaussian, then Pn(E) 
is also Gaussian. If instead these couplings are binary, 
P(Jij) — [S(Jij + 1) + S(Jij — l)]/2, we are lead to a 
"discrete" REM where E is an integer random variable 
with distribution P N (E = -dN + 2k) = 2~ dN B(dN, k). 
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In this expression, B is the usual binomial coefficient giv- 
ing the number of ways to choose k elements out of dN. 
This distribution is roughly Gaussian at large N. The dif- 
ferences between the discrete and continuous models be- 
comes most apparent for the extreme energies and thus in 
the low temperature phase where the lowest levels dom- 
inate the partition function. 

Let us begin with the properties at T = 0. The ground 
state energy Eq grows linearly with N and its variance is 
0(1) in the large N limits . The ground state is unique 
in the Gaussian case whereas it has a strictly positive 
probability of being degenerate in the discrete case. This 
has important consequences. Consider the overlap dis- 
tribution Pj(q) for a given disorder instance; when the 
ground state is non-degenerate, Pj will be a delta func- 
tion at q = 1, while when it is degenerate, Pj will have 
two delta function peaks, one at q = and one at q = 1. 
(By convention, any two distinct levels are taken to have 
zero overlap.) Since the degeneracy does not disappear 
at large N as we show below, P(q), the disorder average 
of Pj, tends towards two delta function peaks, each of 
weight 0(1) as N -> oo. Thus at T = there is RSB in 
the discrete REM but none in the Gaussian REM. 

Now we shall be more quantitative and determine the 
complete behavior of Pj(q) at zero temperature. Indeed, 
if the ground state is g-fold degenerate for a particular 
disorder instance, one has 



Pj(q) = (l-g- 1 )5(q)+g- 1 6(q-l). 



(2) 



The problem then reduces to computing the statistics of 
the integer g. This can be done analytically as follows. 
Denote by Q%(n, g) the probability that the ground state 
energy is —dN + 2n and its degeneracy g. In the discrete 
REM, we have 



Q%(n,g) = B(2 N ,g) p N (n) 



JN 



k=n+l 



(3) 



where we have introduced pjv(n) = Pn{E = —dN + 
2n) to simplify the notation. Now we are interested in 
calculating the probability of the degeneracy 
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n=0 



From that we shall calculate 



Or 1 ) ^O^^-i 
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is negligible (recall that E and n are related by E 
—dN + n). Since 



p N (n - An) = p N (n) 



a(n) 
1 — a(n) 



1 °'ik 



(0) 

holds for any finite integer An, the last factor in eq. (||) 
for a(n) < 0.5 can be estimated as 
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2 -g 



i-j2p N (k) 



fc=0 



2" -g 



[1 - a(n)]p N (n) 



2 -g 



exp 



1 - 2a(n) 
[l-a(n)}[2 N -g]p N (n) 



(a(n) < 0.5). (7) 



1 - 2a(n) 

By substituting this equation into eq. (H), we find 



- JY 



(n,g) 



[ PN (n)2 N ] 9 exp 



[l-a(n)]2 N p N (n) 
1 - 2a(n) 



(8) 



provided that a(n) < 0.5 and g is not of order exp(iV). 

Now let us denote by n* the integer for which the value 
of pfq{n)2 N is the closest to 1. Because eq. (||) is also valid 
for n = n*, pAr(n* + An)2 N for finite An is expressed as 



PN (n* + An)2 iv = C* exp(a*An) \ 1 + O [ — 



where 



a* — log 



1 - a" 



(9) 



(10) 



a* = a(n*), and C* = pN{n*)2 N . Now the condition 
PN(n*)2 N ps 1 determines n* and leads to the following 
equation determining a* in the large N limit: 



(1 - d) log(2) - dlog(l - a) + a* log 



l-a* 



- 0. 
(11) 



The substitution of eq. (S) into eq. (0) leads us to 



which according to eq. (|2|) gives us the disorder averaged 
P(q) of the model. 

Now let us calculate Q%(n,g) for a(n) = 4^ < 0.5. 
We are not interested in the case a(n) > 0.5 because 
the probability that the ground state energy is positive 



- N 



(n* + An, g) « -(C*) 9 exp(#a* An) 



x exp 



[l-a*]C* exp(a*An) 
1 - 2a* 



(12) 
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and these two quantities become equal as N — > oo. 

To obtain a closed form expression for (g~ 1 ), we now 
approximate the summation in eq. (0) by an integral and 
find in the large N limit 

/OO 1 
-oo 9- 



x exp 



[1 - a*]C*exp(a*a;) 



1 - 2a* 



1 / l-2a* \ g 
ga* \ 1 — a* 



(13) 



This formula tells us that the probability that the ground 
state is degenerate is non-zero even in the limit N — > oo 
and that the distribution of g is well behaved, i.e., all 
moments are finite. From eq. (|l3|), we finally obtain 



x \ 1 — x 



a" 



o 



where we have used 



v^ 9 r d *i ( 1 



(N<i)- 



(14) 



(15) 



one associates a random energy. Finally, the energy of a 
state (leaf of the tree) is given by the sum of the energies 
of the L branches connecting it to the root of the tree 
(residing at level 0). 

Just as for the REM, the low temperature phase of the 
GREM is sensitive to the detailed distribution of the en- 
ergies. In the context of our study, we see that the prop- 
erties found for the REM extend to the GREM as follows. 
At T = 0, we need consider only the ground states. At 
each level, there is a strictly positive probability to have 
degenerate lowest energies for the model having discrete 
energies. One thus has a strictly positive probability to 
find any of the possible overlaps when considering ground 
states only. This is to be contrasted with the Gaussian 
model for which P(q) = S(q — 1) when T = 0. Similarly, 
at T > 0, P(q) will be quite sensitive to the nature of the 
energies as long as we stay within the spin glass phase. 
At very low temperatures, the weights of its peaks other 
than at q = 1 will be linear in T in the Gaussian case 
while they will have an exponentially small temperature 
dependence in the discrete case (cf. eq. [l6]). Finally, as 
one approaches the highest of the critical temperatures, 
many branches at each level contribute to the partition 
function and so from there on the detailed distribution 
of energies becomes irrelevant. 



If we apply these expressions to our model when d = 3, 
we find a* = 0.317 and (g' 1 ) = 0.824; the peak in the 
overlap distribution is thus much lower at q = than at 
q — 1. This finishes our analysis at T = 0. 

Consider now what happens at T > 0. The detailed 
nature of P^{E) will remain important as long as the 
partition function is dominated by a finite number of 
energy levels. This happens throughout the whole low 
temperature phase (T < T c ) where the free-energy and 
Pj(q) are not self-averaging. Consider in particular the 
T dependence of P(q) at low temperature. In the Gaus- 
sian case there is a non-zero probability density to have 
a zero energy gap; this leads to a P(q, T) that is linear 
in T as T — > 0. On the contrary, in the discrete case we 
have 



P(q = 0, T) = P(q = 0, T = 0) + 0{e 



(16) 



where AE — 2 is the energy gap of this model. Indeed 
this holds for each disorder instance and so also holds for 
the disorder average. 



III. GENERALIZED RANDOM ENERGY MODEL 

In the GREmI one considers 2 N states that are orga- 
nized in a tree-like fashion, allowing for correlated ener- 
gies. At the first layer of the tree, there are 2^ branches; 
then each such branch gives rise to 2 N2 branches in the 
second layer, etc... The nodes of the last layer are iden- 
tified with states (or spin configurations), and there are 
2 N of them where N = Ni + N 2 . . . + N L . To each branch 



IV. INFINITE-LEVEL GREM 

In a discrete GREM with L-layers, the ground state 
may be degenerate for any finite L, but what happens 
when L — > oo? To study that limit, we set Ni — k, k 
being a fixed integer (say 1 or 2). The model is schemat- 
ically shown in figure |l|. There is then a fixed branching 
factor K = 2 k at each layer, each node generating K 
branches. A random energy e is associated with each 
branch of the tree. The e variables are independent and 
drawn from the same distribution pe(^)- The energy E{i) 
of a state i is given by summing up the e's of the branches 
which lie along the path connecting it to the tree's root O, 
e.g., in figure |], E(i) = ei+£2 + £3 and E(j) = ei + e 2 + e 4 . 
The distance dij of two states i and j is d if their first 
common ancestor arises on the d-th layer counted from 
below, e.g., in figure 0, da = 0, dij = 1 and dik = djk = 2. 
The overlap qij is related to d^ by 



qij = 1 - d^/L. 



(17) 



Thia-jnodel haSrbeen studied in depth by Derrida and 
SpohnllJ (see alsolI3tL3); it can be viewed either as an 
infinite level GREM or as a directed polymer on a dis- 
ordered Cayley tree. When the energy of each branch is 
taken from continuous distribution, the model's thermo- 
dynamics is extremely close to that of the REM. There 
is a critical temperature below which a finite number 
of states dominate the partition function, and this low 
temperature phase exhibits one-step RSB, while the dis- 
tribution P(q) tends towards a delta function at q = 1 as 
T^0. 
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FIG. 1. Construction of the infinite-level GREM with a 
branching factor K — 2. 

Since we are motivated by the question of RSB in the 
3-dimensional ± J EA model, we shall investigate the be- 
havior of P(q) for the Cayley tree model when the ener- 
gies on each branch are discrete. We shall take pe{() to 
consist of a finite number of delta functions. To investi- 
gate the behavior of our model, we shall derive recursion 
formulae for the probability Y(L,d) to find two states at 
a distance less or equal to d by using-the same techniques 
as developed by Derrida and SpohnEl Since q and d are 
related by eq. Jl7|), we can calculate P(q) from Y(L,d). 
By definition, the probability Y(L,d) is given by 



Y(L,d)=^-^ exp[-X(i)-X(i)]. (18) 

ij/dij<d 

In this expression, represents the disorder average, 
X(i) = E(i)/T, Z is the partition function of the sys- 
tem, and L is the number of layers. A standard integral 
representation for Z~ 2 leads to 



/oo 
duF(L,d; u), 
-OO 



(19) 



where 



F(L,d; u) 



= cxp[-2u-c- u Z{L)} e-^W-^O'). (20) 

ij/dij <d 

Collecting terms in the sum that belong to the same sub- 
tree of height d, we have 



ij/dij <d 



-x(i)-xu) 



^exph2A( J B rf )]z(B d ) 2 . (21) 



In this expression, B d is a general branch point in the 
d-th layer counted from below, z(B) is the partition 
function of the sub-tree rooted at a branch point B, 
X(B) ee E{B)/T, and E(B) is the energy of a branch 
point B given by summing up the e's of the branches 



which lie al ong the path connecting B and O. Substitu- 
tion of eq. d2l| ) into cq. ( pp| ) gives 

F(L,d; u) 

= exp[-2u - e-^Z(L)} ^ exp[-2X(B d )}z(B d y. (22) 



From this equation, we find 

F(d,d;u) =H 2 (d;u), 

where 



H n (d;u) ee {e- u z(B d )} n ex V [-e~ u z(B d )}. 



(23) 



(24) 



We can calculate H n (d; u) (including H2 (d; u) appear- 
ing in eq. (23)) by the following recursion formulae. Now 
let us start from the simplest case, i.e., n = 0. By defini- 
tion, 



H (0;u) =exp[-e- u ]. 
Since the z(B d +i) can be expressed as 
z(B d+1 )=Y,o~ X(Bd) z(B d ) 



(25) 



(26) 



and z(B d ) and z(B' d ) are independent if B d 7^ B' d , we 
obtain the recursion formula 



H (d + 1; u) = JJ cxp[-c-«-^(^)z(B d )] 

B d 

= H (d;u) K , (27) 
where we have defined 

g(u) = J dep E (e)g(u + e/T) (28) 

for a general function g{u). The recursion formulae for 
tj ^ are derived by using eqs. ( p5|) and ( f27|) as well as 
the relation 



H n (d;u) = -^-H (d;u). 



(29) 



For example, the recursion formula for H\(d] u) is 
H^d+l-u) 



^-H (d;u) K 
au 



= KH^d-^H^u) 



K-l 



(30) 



Finally, let us derive recursion formulae for F(L,d;u). 
From eq. (p2|), we have 



F(L + l,d;u) 



exp 



-2u - e"" Y e~ x{B ' L) z{B' L ) 
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^ exp[-2X(B d )]z(B d ) 2 
exp [-2u - 2X(B L ) - C - u - x( - B ^ z(B L )] 



e - 2 i x ( B ^- x ( B ^z(B d ) 2 

B d eB L 



x j I exp [— e 



u-X(B' 



>z(B> L 



(31) 



The strategy now is to perform the disorder average in 
two steps: first we average over sub-trees rooted at {Bl} 
when fixing {X(Bl)\; then we average over {X(Bl)}. 
By using eq. (|2|) and eq. (|I|) with n = 0, the disorder 
average in the first step is done as 

F(L + 1, d; u) = ^F(L,d;u + X(B L )) 

B L 
X 



[] H (L;u + X(B' L )). (32) 

B>=£B L 



Note that the only random variables in this equation are 
{X(Bl)}. The disorder average in the next step finally 
leads us to 



F(L + 1, d; u) = KF(L, d; u)H (L; u 



K-l 



(33) 



In summary, the disorder averaged distribution of dis- 
tances Y(L,d) can be computed by the following proce- 
dures: 

i) Calculate iJ 2 (d;w) (=F(d,d;u)) by evaluating nu- 
merically the recursions which are derived by ap- 
plying eq. © to eqs. @ and @. 



ii) Calculate F(L,d; u) by using the recursion eq. (33). 



iii) Compute Y(L,d) by estimating numerically the in- 
tegral in eq. (|l9|). 

For large number of layers we see convergence to the 
large L limit after introducing the (continuous) overlap 
q = 1 — d/L; this gives us the disorder average P(q) or 
equivalently Y(q) = P(q)dq for the infinite tree. 

Just as for the Gaussian case, the discrete model has 
one step RSB below a critical temperature T c . However, 
the energies are discrete and we find that the degeneracy 
of the ground state is maintained as L — > oo and fur- 
thermore its mean saturates to a finite value. Because of 
this, there is RSB at T = 0, not only for < T < T c . 
When executing our recursions, the distribution of over- 
laps rapidly has very narrow peaks at q — and q = 1 
with increasing number of layers, while the probability 
of having overlaps between those two values goes to zero. 
One can thus characterize the mean probability distribu- 
tion of overlaps by the weight of the delta function at 
q = 0. For that, we consider the quantity Y(q = 0.5) 
which gives the weight of the overlaps near q = 0. In 
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FIG. 2. Y(q = 0.5) = / Q P(q)dq as a function of tem- 
perature, for the infinite level GREM (top curve) and for its 
random entropy generalization (bottom curve). (Plots are for 
K — 2 and L = 1000.) 



figure |2j we show the temperature dependence of that 
quantity (top curve). For these data, we used K = 2 and 

p E {e) = 0.255(e) + 0.5£(e - 1) + 0.25<5(e - 2). (34) 

Furthermore, we took L — 1000 which is big enough to 
represent the infinite tree limit at all T except for T close 
to T c . We have checked that the curves near T c converge 
as L increases, and a cusp is created at T c in the limit 
L — > oo, just as the continuous caseEj. But the point we 
want to make here concerns T -C T c : the curve is very 
flat at low temperatures. This is as expected from the 
analogy with the discrete REM (cf. eq. [16]), and is due 
to the presence of a gap above the ground state energy. 



V. A MODEL WITH EXTENSIVE ENTROPY 

In all the models considered so far, a finite number 
of the lowest energy states dominate the partition func- 
tion at sufficiently low temperatures. Thus the entropy is 
0(1) rather than O(N). Clearly, in any physical system, 
the entropy remains extensive as long as T > 0. What 
is the possible importance of such extensivity? Krzakala 
and Martintil argued that extensive entropies will give rise 
to diverging entropic fluctuations, and that these fluctu- 
ations should prevent RSB at T = 0. To test this idea, 
one needs a model with extensive entropies at low tem- 
perature. There are several ways to do this-Ln the context 
of GREM-like models; Cook and—DerridaEj added small 
loops to the tree, while YoshinoEj replaced the random 
energies by random entropies. We keep the tree structure 
but have both energies and entropies. 

We focus again on the infinite tree with a constant 
branching factor K. To each branch we assign a random 
energy e and a random entropy a from distributions Pe{z) 
and ps (c) , respectively; then each leaf has an energy and 
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entropy given by the sum of those terms along the path 
connecting that leaf to the tree's root O. In this model, 
both the energy and the entropy of the states become 
extensive. Now a crucial point is that the model with 
and without entropies can be mapped onto one-another 
when T ^ 0; indeed the two models are identical provided 
the distributions satisfy 



p E (e') = / deda<5(e' - e + Ta)p E {e)p s {a), 



(35) 



where p* E {e') is the distribution of the energy for the 
model without entropy. This condition insures that 
the distribution of the weight assigned to each branch 
(exp(— e/T + a) and exp(— e'/T)) is the same in the both 
cases. Therefore, we can utilize all of the recursions pre- 
viously derived just by changing eq. (f28[) into 



g(u) = J dedap E (e)p s (a)g(u + e/T - a). (36) 

We have investigated the distribution of overlaps in 
this model with entropic fluctuations. We have chosen 
Pe as before and set 



ps(a) =0.5<5(cr) + 0.5<5(cr-2). 



(37) 



We have also chosen the same branching factor as before, 
i.e., K = 2. The inclusion of these entropic fluctuations 
changes the critical temperature, but it also changes the 
qualitative behavio r of RSB. In figure || we show the same 
quantity as before, Y(q — 0.5). We see that the entropic 
fluctuations destroy RSB at T = 0, and P(0) goes to 
zero linearly with T as T — ► 0. One may think from the 
figure that there is no cusp at T = T c in the model with 
entropy. However it can be shown from the mapping 
mentioned in the previous paragraph that the inclusion 
of entropy does not remove the cusp. To check this, we 
have verified numerically that the cusp appea rs slowly as 
L is increased and that the left derivatives of Y(q = 0.5) 
at T = T c tend to a non-zero limit for both models as 
L — ► oo. 

We have also investigated the finite size dependence of 
P(q = 0) at T = 0. From figure ||, we see that this quan- 
tity decreases as L~ 5 . This result can be understood 
from the following simple argument. Imagine a disorder 
instance in which there are only two ground states whose 
overlap is close to 0. The probability for an overlap to 
be small for such a sample is of order 1 if the difference 
of the entropies of the two ground states is of order 1. 
Recall that each ground state has an extensive entropy; 
furthermore, the entropy fluctuations of a ground state 
are necessarily of order L 5 . Then the probability that 
the two states have the same entropy is of order L~ 5 
in the large L limit. From this we conclude that at zero 
temperature, P{q = 0) decreases as L -0 5 as L — > op. 

It is worth pointing out that Krzakala and MartintJ sug- 
gested that in the 3 — d EA ± J spin glass valley-to- valley 
fluctuations in the entropy should grow as L d "/ 2 (d s is the 
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FIG. 3. Y(q = 0.5) at T = vs. L. The parameters for 
these data are K = 2, ps(o') = 0. 55(c) + 0.5<5(<r — 1) and 



eq. ( p^ ) for P e- A function proportional to L 
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fractal dimension of surfaces of large-scale low-energy ex- 
citations) and P(0) should then decrease as L~ ds / 2 . A 
power law decay of P(0) in that model has been found 
numericallyt3'Q The contribution of our study is to show 
that such a phenomenon does indeed occur beyond rea- 
sonable doubt in a particular model. 



VI. DISCUSSION 

Let us compare the model with discrete energies and 
extensive entropies to the one where the disorder vari- 
ables are continuous (be there entropy fluctuations or 
not). In both cases, P(q = 0) goes to zero linearly with 
T as T — > because the lowest free- energy states are 
non-degenerate and there is no gap. We may then sum- 
marize what we have found by saying that state-to-state 
entropy fluctuations of order y/L effectively remove both 
degeneracies and gaps. Of course at T = 0, the degener- 
acy in the energy is important, but our point is that the 
free-energies are still not degenerate: the contribution 
of different ground states to the partition function are 
wildly different simply because their entropy differences 
diverge as \f~L. This is the mechanism behind no RSB at 
zero temperature. If we consider now T > 0, we note that 
the energies of the states dominating the partition func- 
tion are far above the lowest energy; the system selects 
the states with the minimum free-energies, F = E — TS; 
even though both E and S are integers, generically F is 
not and so any fine structure of E and S is washed out. 

Our work has focused on the infinite tree model be- 
cause it is tractable, but we expect the conclusions to 
be quite general when there is one-step RSB at T > 0. 
An open question concerns of course the case of systems 
having continuous RSB. One way to address this is to 
generalize our model so that it exhibits continuous RSBL 
this can be done just as for the Derrida-Spohn modelcHl 
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by making the distributions of energies depend on the 
level of the tree. With such a modification, the model 
having no entropy but discrete energies will have con- 
tinuous RSB at T = 0, whereas if random entropies are 
assigned to each branch the entropic fluctuations will re- 
store replica symmetry at T = 0. 

One may also ask what effect would Jij = ±1 have 
on the Sherrington-KirkpatrickEZI (SK) model, i.e., would 
such discreteness give rise to RSB at T = 0, in contrast to 
what happens in the Gaussian case? The crucial prop- 
erty that can make RSB possible at T = is ground 
state degeneracy; the gap in such a discrete SK model 
of N spins is l/\J~N if the J^s are rescaled so that the 
model has a thermodynamic limit. The valley-to- valley 
energy differences being O(l), the probability of an ex- 
act ground state degeneracy should be 0(1 /\N), so we 
do not expect RSB at T = here. In fact, it seems un- 
likely that the discrete and continuous SK models have 
any differences in the thermodynamic limit because the 
local fields on the spins effectively become continuous in 
the large N limit. 

Finally, let us note that the most controversial case of 
zero-temperature RSB arises in the 3 — d EA ± J spin 
glass; its breaking of replica symmetry at T > is prob- 
ably continuous, though even that is subject to debate. 
It would thus be very useful to puxsuc. this issue further 
to resolve the contradictory resultsE^ruCm obtained so far. 
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